Atomic-scale relaxation dynamics and aging in a metallic glass probed by X-ray 

photon correlation spectroscopy 
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We use X-Ray Photon Correlation Spectroscopy to investigate the structural relaxation process 
in a metallic glass on the atomic length scale. We report evidence for a dynamical crossover between 
the supercooled liquid phase and the metastable glassy state, suggesting different origins of the relax- 
ation process across the transition. Furthermore, using different cooling rates we observe a complex 
hierarchy of dynamic processes characterized by distinct aging regimes. Strong analogies with the 
aging dynamics of soft glassy materials, such as gels and concentrated colloidal suspensions, point 
at stress relaxation as a universal mechanism driving the relaxation dynamics of out-of-equilibrium 
systems. 
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Glasses are usually defined as liquids that are trapped 
in a metastable state from which they slowly evolve 
toward the corresponding equilibrium phase 0, Q. 
Although aging is known since centuries, a clear picture 
of the dynamics in the glassy state is still missing 0, 0] . 
Most of the experimental information available on aging 
concerns macroscopic quantities, such as viscosity or 
elastic moduli jH-ll|, or focuses on dielectric relaxation 
12H171] , a quantity that is often difficult to relate directly 
to the particle-level dynamics. From these measure- 
ments, a characteristic time for the evolution towards 
equilibrium can be extracted, but no direct information 
on the connection between aging and the underlying 
microscopic dynamics is available. By contrast, a full 
understanding of aging requires a detailed description 
of the particle-level dynamics, and in particular of the 
structural relaxation time r, the characteristic time 
for a system to rearrange its structure on the length 
scale of its constituents. While the structural relaxation 
has been widely investigated in the liquid phase, only 
a few studies report on the beh aviour of r below the 
glass transition temperature, T g 
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Molecular 

dynamics simulations show that t increases linearly or 
sub- linearly with the waiting time, t w , and suggest the 
possibility of rescaling the measured quantities on a 
single master curve (time-waiting time superposition 
principle) fl8j . Studies in this direction, however, have 
led to debated results d E H, E3, 113 ■ Consequently, 
several key questions remain unanswered: what is the 
fate of the structural relaxation process when the system 



falls out of equilibrium in the glassy state? How does 
it depend on the thermal history and the waiting time? 
What physical mechanism is responsible for structural 
relaxation? 

Here, we address these questions by presenting an 
experimental investigation of the structural relaxation 
on the atomic length scale in a metallic glass former, 
in both the supercooled liquid and the glassy state, as 
a function of temperature and waiting time, and for 
different thermal histories. We use X-Ray Photon Cor- 
relation Spectroscopy (XPCS) to study Mg 6 5Cu25Yi , a 
well-known glass former with a relatively low T g ~ 405 
K and a stable supercooled region. 

The experiments were performed at the bcamlinc ID10 
of the European Synchrotron Radiation Facility. A 
partially coherent beam with an energy of 8 keV was 
selected by rollerblade slits opened to 10 x 10 fim, placed 
~0.18 m upstream of the sample. Mg6sCu25Yio ribbons 
were placed in a furnace, providing a temperature 
stability of 0.1 K. The scattering intensity was collected 
in transmission geometry by an IkonM charge coupled 
device from Andor Technology (1024 x 1024 pixels, 
13 x 13 \ivn pixel size) placed 0.68 m from the sample. 
Sets of ~ 3000 images were taken with 5 s or 7 s 
exposure time per frame and were analyzed following the 
procedure described in Ref. [2(| . The quantity accessible 
in a XPCS experiment is the intensity auto-correlation 
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denotes the ensemble average over all the pixels of the 
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detector and (...) is the temporal average. g2(q,t) can 
be related to the intermediate scattering function f(q,t) 
through the Siegert relation g%(q, t) = 1 + B(q) \f(q, t)\ , 
with q the scattering vector and B{q) a setup-dependent 



parameter 2l|. By measuring the signal of all pixels, the 



intensity distribution is indeed sampled correctly allow- 
ing the extension of the Siegert relation to non-ergodic 
systems [22 - 24 1 . We perform XPCS measurements 
at a wave-vector qo = 2.56 A" 1 corresponding to the 
maximum of the static structure factor, thereby probing 
directly the dynamics at the inter-particle distance 
2n/q ~ 2 A. 

Our sample was produced by melt spinning with 
an extremely fast quench (10 6 K/s) from the high- 
temperature melt down to T = 300 K (see SI). We first 
measure this "as-quenched" glass (AQ-g) by performing 
several isothermal runs at increasing values of T, up 
to T g . We then explore the supercooled liquid region 
(T > T g ), before slowly cooling again the system below 
T g ("slowly-cooled glass", SC-g). A second as-quenched 
glass (AQ-g2) is also measured, directly heated to 415 
K without performing intermediate steps. Between 
isothermal runs T was changed at a rate of 1 K/min. 
Figure Q] shows a selection of normalized g2 (t) measured 
for temperatures above and below the calorimetric T g for 
the SC-g. The data are reported together with the best 
fit obtained by using the Kohlrausch- Williams- Watt 
expression g2{t) = 1 + Aexp(— 2(i/r)' 3 ), with r the 
structural relaxation time, ft the shape parameter 
and A = A(q,T) = B(q) * f q (T), where f q (T) is the 
non-ergodicity factor appearing in the expression of 
the final decay of the intermediate scattering func- 
tion, f(q,t) = f q (T) exp(— t/r)P. A first important 
observation is that the dynamics is not frozen in the 
glassy state: gi (t) — 1 displays a full decay to zero, 
whose characteristic time r increases with decreasing 
temperature. Thus, structural relaxation occurs even if 
the system has not fully equilibrated A full decay is 
not completely reached at lower temperatures due to the 
large relaxation time of the sample. As discussed below, 
in the case of AQ-g the dynamics is much faster and we 
do observe a full decay at all T . More surprising, we find 
a remarkable equilibrium/out-of-equilibrium dynamical 
crossover at T g as signalled by a sharp change of the 
shape of g2(t), from stretched (ft = 0.88 < 1) in the 
supercooled liquid to "compressed" (ft ~ 1.3 > 1) in the 
glass. As it will be discussed later, the curves in the 
glass depend on waiting time and thermal history, which 
is not the case for data collected above 408 K. 
Stretched exponential functions are a universal feature 
of supercooled liquids 0, reflecting a wide distribution 
of relaxation times, generally attributed to the existence 
of dynamical heterogeneities 




25J, |26|. Differently, the 



compressed behaviour found in our structural glass is 
totally unexpected and cannot be explained as resulting 
from a distribution of single exponential relaxation 
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FIG. 1: Temperature dependence of correlation functions 
measured for qo = 2.56 A -1 by means of XPCS. Lines are 
the best fits using a KWW model function. The black arrows 
emphasize the dynamical crossover from stretched (ft < 1) 
to compressed (ft > 1) exponential functions below 408 K. 
The curves in the glassy state correspond to the SC-g and are 
measured on lowering T from the liquid phase, b, Tempera- 
ture dependence of the corresponding shape parameter. The 
line indicates the T g obtained from calorimetric measurements 
(see SI). 

times, as it would be within the heterogeneous dynamics 
scenario. Indeed, this scenario could explain a growth 
of ft in the out-of-equilibrium state, but only up to 

ft < 1 [H, |27H3(j. The idea would be that below T a 
only the regions characterized by a faster-than-average 
dynamics contribute to the decay of the correlation 
function, leading to a narrower distribution of relaxation 
times and thus to a ft closer to, but still smaller than, 
,5 = 1, the limiting case of a single exponential decay. 
An alternative explanation for the growth of ft might be 
that the relaxation time evolves during the measurement 
of c?2 (*) , altering the shape. We rule out this hypothesis 
by verifying that ft does not depend on the time interval 
over which <?2 (t) is averaged, and by noting that the 
slowing down of the dynamics below T g would lead to 
decrease the measured value of ft, rather than increasing 
it 0. 

Compressed correlation functions have been observed 
in out-of-equilibrium soft materials, such as concen- 
trated colloidal suspensions and polymeric gels, and for 
nanoparticle probes in glass former matrix 3lM38|. In 
the latter case a similar stretched/compressed crossover 
was reported; however, those experiments probed length 
scales much larger than that of the glass former units 
and located the crossover in the liquid phase, well 
" 38 



above T 37 



By contrast, our experiments provide 
strong evidence of compressed exponential relaxations 
as soon as the system falls out of equilibrium below T g . 
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An increase of (3 for T < T g was also reported in the 
investigation of the dielectric spectra of some polymeric 
glasses ljj However, in these systems /? increases by 
not more than ^10% with respect to the corresponding 
liquid value, and always remains less than one. In our 
case, not only does /? increase by more than ^40% in 
the SC-g (Figure [Tp), but it is also nearly temperature 
independent, in contrast with the prediction of a recent 
theoretical work pjoj ]. 

Since diffusive processes driven by thermal energy lead 
to an exponential (or stretched exponential) structural 
relaxation, the decay observed here points to a different 
origin of the microscopic dynamics, unaccounted for in 
previous studies 

HUM- 

We propose internal stress 
relaxation as the main physical mechanism responsible 
for the dynamics in the glass. Indeed, in experiments 
on soft materials a similar compressed exponential 
relaxation has been associated with ballistic rather 
than diffusive motion due to internal stresses stored in 
the system [3, These stresses can be described 

by progressive irreversible rearrangements that act as 
local sources of dipolar stress leading to a compressed 
exponential shape of the dynamic structure factor (35j . 
It is well-known that the fast quenching used to produce 
metallic glasses leads to internal stresses or excess frozen 
in free volume 0, I39H42T ]. which may be released by 
heating the system (see Dynamic Mechanical Analyzer 
(DMA) and Differential Scanning Calorimetry (DSC) 
data in SI). 

The temperature dependence of the structural relaxation 
times is shown in Figure [2] Data corresponding to three 
different thermal paths are presented (full symbols): i) 
diamonds are taken from the curves reported in Figure 
[1] corresponding to the supercooled liquid and to the 
SC-g; ii) squares correspond to the AQ-g; hi) the star 
corresponds to AQ-g2. Macroscopic measurements 
obtained by viscosimetry, DMA and DSC are also 
shown as empty symbols. The viscosity data 43j have 



been rescaled to the relaxation time values measured 
at higher temperatures with DMA. Data calculated 
from DSC measurements are also reported (see SI). The 
XPCS data taken above 403 K follow the macroscopic 
equilibrium liquid line and correspond to the stretched 
exponential correlation curves reported in Figure [1] 
typical of supercooled liquids. All the points below 
the equilibrium line are instead in the glassy state. 
They come from compressed exponential curves and 
display aging. In the glassy state t shows a strong 
cooling-rate-dependent departure from the behaviour 
observed at equilibrium, being the dynamics much faster 
in the AQ-g than in the SC-g. 

On increasing T in the AQ-g, the dynamics gets 
smoothly faster with T, until ~ 393 K where r in- 
creases sharply with increasing temperature, before 
decreasing again above 408 K, once the equilibrium 
curve is reached. This surprising increase is likely just 
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FIG. 2: Structural relaxation times measured with XPCS. 
The dotted-dashed line indicates the T g for the SC-g. We 
report data for each glass at the same waiting time t w ~ 6000 
s for the AQ-g and t w ~ 9000 s for the SC-g after temperature 
equilibration. The filled circles correspond to the AQ-g at the 
same t w ~9000 s as for the SC-g, for a consistent comparison. 



a consequence of annealing during isothermal runs. 
Here the system slowly relaxes its internal stresses 
during heating and loses memory of its original thermal 
treatment, thus falling on the upper equilibrium curve 
6, 44]. By contrast, if the system is continuously heated 
up to a higher temperature without having the time 
to rearrange (AQ-g2, in Figure [2]), this feature is not 
observed, and the glass meets its corresponding liquid 
at higher temperatures. 

When the system falls out of equilibrium it starts 
aging. The age dependence of the dynamics is directly 
captured by XPCS through the analysis of the two-time 
correlation function Cfafa) (Figure [3]). This quantity 
corresponds to the instantaneous correlation between 
two configurations at times t\ and t-i- Lines parallel 
to the diagonal from the lower left to the upper right 
corners are iso-lag time lines (t = ii-£2=constant), while 
the sample a ge v aries along such lines according to t age 



(^i+^)/2 [2l|. The slowing down of the dynamics 



is illustrated by the broadening of the yellow-red ridge 
along t=0 as the age increases. Since the evolution of the 
dynamics due to aging is slow, we time-average C(t\,t2) 
over short time windows to obtain age dependent 
correlation functions g2(t,t w ). Here t w = £o+(t/-t»)/2, 
to being the delay of the start of the measurement from 
temperature equilibration, and tf and ti the final and 
starting times of the interval considered for averaging 
Cfafy). Examples of different g2(t,t w ) measured for 
different ages are reported in Figure [3j3. On increasing 
t w , the dynamics shifts towards longer time scales. 

Figure [Ik. reports the t w dependence of (3 in the 
AQ-g for T >318 K. We find a temperature- and 
t^-independent value /3=1.5, which is ^70% larger than 
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FIG. 3: a, Two time correlation function measured in the 
AQ-g at T=318 K by means of XPCS. b, Corresponding 
92{t,t m ) calculated for different sample ages. From left to 
right t u) =4680 s, U=9840 s and £^=12960 s. The data are 
reported together with the best fit curves with a KWW model 
function. 

in the supercooled liquid. This value is also larger 
than that for the SC-g (j3 ~ 1.3), confirming that 
deviations from diffusive behaviour are enhanced by a 
faster quench, consistent with the internal stress picture. 
The value /3=1.5 has been often observed in studies of 
soft materials 3l|, [33J, |3J] and it is also in agreement 
with the proposed models [34, IH. Only at 403 K does 
/3 decrease at very large t w , presumably indicating that 
the sample gets closer to equilibrium. 

Surprisingly, the aging differs from the linear or 
sublinear laws reported in the literature and exhibits 
remarkable scaling properties. For each temperature r 
increases with t w exponentially, and the growth rate 
does not depend on T. This suggest that sets of data 
taken at the same temperature may be modelled by the 
empirical law 



t(T, t w ) = T (T)exp(t w /T*) 



(1) 



and rescaled onto a master curve by plotting t/tq vs. t w , 
as shown in Figure 0}}, where i*=7940±70 s independent 
of T. An exponential growth of r with sample age 
has bee n p reviously reported in a few soft materials 
4^-47|: our results show that this unexpected 
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behavior may be more general than previously thought. 
The constant shape parameter and the perfect scaling 
of T(T,t w )/T (T) lead moreover to the validity of a 
time-waiting time-temperature superposition (TTTS) 
principle in the AQ-g and therefore to the possibility 
of obtaining a single master curve plotting 52 (t) as a 
function of the reduced time tj '(T(T,t w ) /t (T)). 
The exponential aging persists even close to T g , when j3 
starts to decrease with increasing t w . Here the system 
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FIG. 4: a, Sample age and temperature dependence of the 
shape parameter /3 measured with XPCS in the AQ-g. b, 
Corresponding structural relaxation times normalized by their 
initial value to(T). Data corresponding to temperatures where 
only one curve has been measured are reported only in panel 
a (empty symbols), c, Structural relaxation times measured 
at T = 403 K for both AQ-g (circles) and AQ-g2 (stars) nor- 
malized by their initial value to and reported as a function 
of t w rescaled by the parameter t* obtained from the analy- 
sis of the fast aging regime. The data are reported together 
with the corresponding equilibrium values (red and green ar- 
row, respectively) obtained by rescaling for to the equilibrium 
value reported in Figure [2] 

enters an intermediate aging regime where both r and 
vary with age. This is the case for the AQ-g data 
taken at 403 K (Figure 0^) and for all data taken in the 
SC-g (see SI). Generally, one expects the exponential 
aging to stop at one point when the materials is par- 
tially annealed. To test this hypothesis, we lower the 
temperature of the AQ-g2 glass to 403 K, after partially 
annealing it at 415 K without equilibration (Figure 
0). At 403 K, the AQ-g2 data display the same fast 
aging as for the AQ-g, although with a higher starting 
value To,AQ-g2(T)=390 ± 20 s and a slightly different 
i*AQ-g2=5600±200 s, due to the different thermal 
history. Differently from the AQ-g, upon increasing 
t w the AQ-g2 glass enters a second dynamical regime, 
where aging slows down so much that it becomes 
undetectable within our experimental window. Here r 
is most likely reaching equilibrium sublinearly, as in the 
case of polymeric glasses and colloids 0, 0, Q . Note 
that here j3 ~ 1.15 is t w independent and still ~ 30% 
higher than in the corresponding liquid, confirming that 
the sample is not fully equilibrated. 
In conclusion, we have scrutinized the dynamics of 
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a metallic glass on the atomic length scale and for 
different cooling rates and thermal paths. We find a 
universal dynamical crossover between the supercooled 
liquid phase and the glassy state where the shape of the 
correlation function changes sharply when crossing T g . 
We have proposed that internal stress relaxation is the 
mechanism driving the dynamics in the glass. Consis- 
tently, we find that deviations from the diffusion-driven 
dynamics typically observed in the supercooled state are 
more important for samples where internal stresses are 
likely to be larger, e.g. samples submitted to a faster 
quench and at earlier ages. As the sample ages, internal 
stress are presumably relaxed, but in a way that depends 
on thermal history, leading to distinct aging regimes. 
Different aging regimes and compressed correlation 
functions have been reported in the past for gels and 
colloids . thus suggesting a universal phenomenology 

umiii-ii^. 
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